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Abstract. Photospheric parameters and chemical composition are determined 
for the single-hned chromosphericahy active RS CVn-type star 29 Draconis 
(HD 160538). From the high resolution spectra obtained on the Nordic Optical 
Telescope, abundances of 22 chemical elements, including the key elements such 
as ^^C, ^'^C, N and O, were investigated. The differential line analysis with the 
MARCS model atmospheres gives Tcs = 4720 K, log g = 2.5, [Fe/H] = -0.20, 
[C/Fe] = -0.14, [N/Fe] = 0.08, [O/Fc] = -0.04, C/N = 2.40, ^^C/^^C = 16. 
The low value of the ■^'^C/^^C ratio gives a hint that extra mixing processes 
in low-mass chromosphericahy active stars may start earlier than the theory of 
stellar evolution predicts. 

Key words: stars: RS CVn binaries, abundances - stars: individual (29 
Draconis = HD 160538)) 

1. INTRODUCTION 

In this work we continue a detailed study of photospheric abundances in RS 
CVn stars (Tautvaisiene et al. 2010, hereafter Paper I) with the aim to get new 
data on chemical evolution of the photospheric and coronal abundance patterns 
and mixing processes in chromosphericahy active stars. Our aim is to investigate 
abundances of more than 20 chemical elements, including ^^C, ^'^C, N, O and 
other mixing sensitive species. We plan to investigate correlations between the 
abundance alterations of chemical elements in star atmospheres and their physical 
macro parameters, such as the speed of rotation and the magnetic field. 

Here we present the results of the analysis for the RS CVn binary 29 Draconis 
(HD 160538) consisting of a KO HI star {V « 6.64 mag) and a white dwarf. Since 
the primary dominates the spectrum, the photospheric abundances can be rela- 
tively well determined, and thus the system is a good candidate for studying the 
coronal-to-photospheric abundance patterns. 

Hall et al. (1982) were the first to suggest that 29 Draconis is a RS CVn binary. 
Its late spectral type (KOIII according to Bidelman 1954 or K2III according to 
Roman 1955), the Can H & K emission in its spectrum (Bidelman 1954), radial 
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velocity values by Abt & Biggs (1972) and photometric variability reaching = 
0.12 mag, as well as a rotation period of 31.5 days found by Hall et al., allowed 
the authors to conclude that 29 Dra is a RS CVn binary. 

Ultraviolet spectra of 29 Dra were observed by Fekcl & Simon (1985). It 
was found that the absorption lines are slightly broadened by stellar rotation 
with ?;sinz = 8 ± 2 kms~^. The Ha absorption line was weak, as if filled in 
by emission, and variable. In five observations obtained over a period of foiu 
months, the equivalent width ranged from 0.1 to 0.8 A. Mgii h and k lines were 
in emission. The study by Fekel & Simon has revealed a white dwarf companion 
of 29 Dra, whose ultraviolet spectrum is matched best to the energy distribution 
of a Teff = 30 000 K and log 5 = 8 model atmosphere. The mass and radius of the 
white dwarf were found to be 0.55 A'/q and 0.012 i?Q, respectively. A radius of 
the primary component was evaluated to be more than 5 Rq. Fekel et al. (1993) 
derived t;sini = 7.2 kms~^, the orbital period 903.8 days and the parameters for 
the companion star. They suspected that 29 Dra can be a barium-like star. 

However, high-resolution spectral analyses by Berdyugina (1994) and Zacs et 
al. (1997) have ruled out the; hypothesis that 29 Dra is a barium-like star. 

In further investigations of 29 Dra, Fekel (1997) determined a more accurate 
value of the radius of the primary component, 4.2 Rq, and the projected rotational 
velocity vsini = 6.7 kms^-'^, which was used in our work. 

A long-term photometry of 29 Dra, which lasted about 28 years, was recently 
reported by Zboril & Messina (2009). They found that this star has three magnetic 
cycles (20.3, 11.1 and 7.6 yr). Modeling of the light curve indicated that the spots 
cover from 4 to 10% of total stellar surface. The photometric period was found to 
be variable in the range 26.4-31.8 d, suggesting differential rotation. 

In our recent investigation of A And (Paper I), we found a low ratio of ^^C/^^C 
= 14, which is lower than is predicted by the first dredge-up theory for a first- 
ascent giant with the parameters of A And. Masses and luminosities of 29 Dra 
and A And are very similar. In this work, it was interesting to check a hint that 
extra mixing processes in low-mass chromospherically active fast rotating stars 
may start earlier than in non-active stars, i.e., below the bump of the luminosity 
function of red giant branch. 

2. OBSERVATIONS AND THE METHOD OF ANALYSIS 

The spectra of 29 Dra were obtained in 1999 August on the 2.56 m Nordic 
Optical Telescope using the SOFIN echelle spectrograph with the optical camera, 
which provided a spectral resolving power of i? « 80 000, for 26 spectral orders, 
each of ~ 40 A, in the spectral region from 5000 to 8300 A. The details of spectral 
reductions were presented in Paper I. 

In the spectra we selected 148 atomic lines for the measurement of equivalent 
widths (Table 1) and 17 lines for the comparison with synthetic spectra. The 
spectra were analyzed using a differential model atmosphere technique described 
in Paper I. Here we present only some details. 

2.1. Atmospheric parameters 

Initially, the effective temperature, Teff, of 29 Dra was derived and averaged 
from the intrinsic color indices {B — V)q and {b—y)o using the corrected calibrations 
by Alonso et al. (1999). The values of color indices B — V = 1.043 and b — y = 
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Table 1. Measured equivalent widths of lines, EW, in the spectrum of 29 Dra. 
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0.664 were taken from van Leeuwen et al. (2007) and Hauck & Mermilliod (1998), 

respectively. A small dcreddcning correction of Eb-v = 0.02, estimated using the 
Hakkila et al. (1997) software, was taken into account. 

The agreement between the temperatures deduced from the two color indices 
was quite good, with the difference 70 K only. No obvious trend of the Fe I abun- 
dances with the excitation potential was found (Figure 1). 

The surface gravity log 5 was found by adjusting the model gravity to yield the 
same iron abundance from the Fei and Fe 11 lines. The microturbulent velocity Vt 
value corresponding to a minimal line-to-line Fe I abundance scattering was chosen 
as a correct value. Consequently, [Fe/H] values do not depend on the equivalent 
widths of lines (Figure 2). 




Fig. 1. The [Fei/H] abundance values versus the lower excitation potential xi. 
The mean abundance ([Fei/H]= —0.20 dex) is shown as a dotted line. 
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Fig. 2. The [Fci/H] abundance values versus the equivalent widths. The mean 
abundance ([Fei/H]= —0.20 dex) is shown as a dotted line. 



2.2. Mass determination 

The mass of 29 Dra was evaluated from its effective temperature, luminosity 
and the isochrones from Girardi et al. (2000). The luminosity log(L/L0) = 1.44 
was calculated from the Hipparcos parallax n = 9.66 mas (van Leeuwen 2007) and 
V = 6.64 mag, the bolometric correction calculated according to Alonso et al. 
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(1999), and the above mentioned Eb-v 
was found. 



0.02. The mass of 29 Dra ~l.2M, 
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2.3. Spectrum syntheses 

The method of synthetic spectra was used to determine carbon abundance 
from the C2 Hne at 5135.5 A. The interval 7980-8130 A, containing strong i^C^^N 
and -'^^C-'^^N features, was used for the nitrogen abundance and ^^C/^^C ratio 
determinations. The ^^C/^^C ratio was determined from the (2,0) ^'^C^^N feature 
at 8004.7 A. All log gf values were calibrated to fit to the solar spectrum by Kurucz 
(2005) with solar abundances from Grevesse & Sauval (2000). In Figure 3 we show 
several examples of synthetic spectra in the vicinity of ^^C^^N lines. 
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X(A) 



Fig. 3. Synthetic spectrum fits to ^^CN lines. The observed spectrum is shown as 
a solid line. The dashed lines are synthetic spectra with [N/Fe] = 0.02, 0.12 and 0.22 dex. 

The oxygen abundance was determined from the forbidden [O i] line at 6300.31 A 
(Figure 4) with the oscillator strength values for ^^Ni and ^°Ni from Johansson et 
al. (2003) and the values log gf — —9.917 obtained by fitting to the solar spectrum 
(Kurucz 2005) and log Aq = 8.83 taken from Grevesse & Sauval (2000). 

The abundance of Nai was estimated using the line 5148.84 A which due to 
rotational broadening is blended by Nil line at 5148.66 A. These two lines are 
distinct in the Sun, so we were able to calibrate their log gf values using the 
solar spectrum. However, the sodium abundance value in our study is affected by 
uncertainty of the nickel abundance determination, originating from the Equivalent 
Widths method. Fortunately, the line-to- line scatter of [Ni/H] determination from 
24 lines of Nil was as small as 0.04 dex. 

For the evaluation of Zri abundance the lines at 5385.13 A, 6127.48 A and 
6134.57 A were used. Evaluation of Yii abundance (Figure 5) was based on 
5402.78 A, Prii on 5259.72 A La II on 6390.48 A, Cell on 5274.22 A and 6043.38 A 
and Ndii on 5276.86 A lines. 

The abundance of Euii was determined from the 6645.10 A line (Figure 6). The 
hyperfine structure of Euii was taken into account when calculating the synthetic 
spectrum. The wavelength, excitation energy and total log gf = 0.12 were taken 
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6300.2 
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6301.0 
X{k) 

Fig. 4. Synthetic spectrum fit to forbidden [O i] line at 6300 A. The observed 
spectrum is shown as a solid line. The dashed lines are synthetic spectra with [0/Fe] = 
-0.14, -0.04 and 0.06 dex. 

from Lawler et al. (2001), the isotopic meteoritic fractions of ^"''^Eu. 47.77%, and 
-'^■''■^Eu, 52.23% and isotopic shifts were taken from Biehl (1976). 

Due to the rotation of RS CVn stars, their spectral lines are broadened, so it is 
important to use a correct value of v sin i in the synthetic spectrum calculations. 



Wc used f sini = 6.7 kms 
spectral lines of 29 Dra. 



from Fekel (1997), which fits well the profiles of 



2.4- Estimation of uncertainties 

The sources of uncertainty were described in detail in Paper I. The sensitivity 
of the abundance estimates to changes in the atmospheric parameters for the 
assumed errors (± 100 K for T^s, ±0.3 dex for log g and ±0.3 km s~^ for vt) is 
illustrated in Table 2. It is seen that possible parameter errors do not affect the 
abundances seriously; the element-to-iron ratios, which we use in our discussion, 
are even less sensitive. 

The scatter of the deduced line abundances a, presented in Table 3, gives 
an estimate of the uncertainty due to the random errors, e.g., in the continuum 
placement and the line parameters (the mean value of a is 0.05 dex). Thus the 
uncertainties in the derived abundances originating from the random errors are 
close to this value. 
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Fig. 5. Synthetic spectrum fit to the Ndii fine at 5276.806 A, Yii at 5402.78 A 
and La II at 6390.48 A. The observed spectrum is shown as a soHd hne. The dashed fines 
are synthetic spectra with [Nd/Fe] = 0.08, 0.18 and 0.28 dex, [Y/Fe] = 0.01, 0.11 and 
0.21 dex and [La/Fe] = 0.07, 0.17 and 0.27 dex, respectively for Ndii, Yii and Laii fines. 
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Fig. 6. Synthetic spectrum fit to the Eu ii fine at 6645.10 AThe observed spectrum 
is shown as a sofid fine. The dashed fines are synthetic spectra with [Eu/Fe] = -0.03, 
0.07 and 0.17 dex. 



Since the abundances of C, N and O are bound together by the molecular 
equilibrium, wc have also investigated how an error in one of them typically affects 
the abundance determination of others. 

A[0/H] = 0.10 causes A[C/H] = 0.04 and A[N/H] = 0.03; A[C/H] = 0.10 
causes A[N/H] = -0.12 and A[0/H] = 0.02. A[N/H] = 0.10 has no effect on both 
carbon and oxygen abundances. 
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Table 2. Sensitivity of the abundances to changes of the atmospheric parameters. 
The table entries show the effects on the logarithmic abundance relative to hydrogen, 
A [A/H]. 
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Table 3. Element abundances relative to hydrogen, [A/H]. a is a standard 
deviation in the mean value due to the line-to-line scatter within the species. 
N is the number of lines used for the abundance determination. 
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3. RESULTS AND DISCUSSION 

As a result, for 29 Dra we found the following atmospheric parameters: Teg = 
4720 K, \ogg = 2.5, -^t = 1-4 kms"!, [Fe/H] = -0.20, [C/Fe] = -0.14, [N/Fe] = 
0.08, [0/Fe] = -0.04, as well as the ratios C/N = 2.40 and ^'^C/'^^C = 16. The 
element abundances [A/H] and the a values (the line-to-line scatter) are listed in 
Table 3 and displayed in Figure 7. 

The effective temperature and surface gravity values for 29 Dra given in dif- 
ferent papers are quite similar. Tcff = 4700 K was determined by Randich et al. 
(1994), 4800 K by Berdyugina (1994) and Zacs et al. (1997) and 4720 K in the 
present paper. Our value of log g = 2.5 is in agreement with the result of Zacs et 
al. (1997). Randich et al. (1994) have found the value larger by 0.1 dex. From the 
wings of MgH lines, Berdyugina (1994) has found log g = 2.7. 

The [Fc/H] value of —0.20 dex determined in the present paper is close to the 
value of Zacs et al., —0.24 dex. A slightly higher value, —0.05 dex, was found by 
Berdyugina and an exceptionally low value of —0.8 dex by Randich et al. 
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In Figure 7, the element to iron ratios for 29 Dra are compared with the 

results of other investigations. The clement to iron ratios for the iron peak, s- and 
r-process elements are solar. Abundances of a-elements are enhanced by about 
0.1 dex. 
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Fig. 7. Element abundances for 29 Dra, as determined in this work (filled dia- 
monds), by Zacs et al. (1997, triangles), Berdyugina (1994, circles) and Randich et al. 
(1994, asterisks). 



3. 1 . The ratios C/N and^'^C/^^C 

The evolutionary sequences in the luminosity versus effective temperature di- 
agram by Girardi et al. (2000) show that 29 Dra with its luminosity log(L/_L0) = 
1.44 is a first ascent giant lying slightly below the red giant sequence bump at 
log(_L/LQ) = 1.6 (Charbonncl & Lagarde 2010). According to the mentioned 
model of mixing, carbon and nitrogen abundances in 29 Dra should be altered 
only by the first dredge-up. The ratio ^^C/^'^C for the first ascent giants with the 
mass of 29 Dra (1.2 M©) should be lowered to the value of about 26 (Charbonnel 
& Lagarde 2010). However, the ratio found in this paper is 16. 

In Figure 8, we compare C/N and ^'^C/^^C ratios of 29 Dra with the theoret- 
ical models including extra mixing. The model, called 'cool bottom processing' 
(CBP), was proposed by Boothroyd & Sackmann (1999) and the model, called 
'thcrmohalinc mixing' (TH), was proposed by Charbonncl & Lagarde (2010). The 
position of 29 Dra in Figure 8 indicates that its ratio of carbon isotopes has been 
altered by extra mixing. The low value of the ^^C/^^C ratio was also found in 
A And (Paper I). 
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Fig. 8. Comparisons of C/N and ^^C/^'^C ratios in 29 Dra (filled diamonds) and 
A And (open diamonds, Paper 1) with the theoretical predictions explained in the text. 



The abundance of lithium is also very sensitive to mixing. During the first 
dredge- up, for a star of the mass of 29 Dra, Li abundance drops to logA(Li) = 
1.45 (Charbonncl & Lagarde 2010). However, the available determinations of Li 
abundance in 29 Dra show lower abundances. Berdyugina (1994) finds logyl(Li) = 
0.6 and Randich et al. (1994) find logA(Li) < 0. The low abundance of Li is 
compatible with the thermohaline extra-mixing model by Charbonnel & Lagarde 
(2010), which predicts the lowering of Li abundance to the value of about 0.07 dex. 
The observed Li abundances in A And are even lower (Savanov & Berdyugina 1994; 
Randich et al. 1994; Mallik 1998). 

Thus, both A And and 29 Dra, investigated by now in our program, give a hint 
that extra-mixing processes may start acting in these low-mass chromospherically 
active fast rotating stars slightly earlier than at the bump of the red giant sequence 
in non-active stars. 
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